Variational Analysis for Photonic Molecules 
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A new type of artificial molecule is proposed, which consists of coupled defect atoms in photonic 
crystals, named as photonic molecule. Within the major band gap, the photonic molecule confines 
the resonant modes that are closely analogous to the ground states of molecular orbitals. By 
employing the variational theory, the constraint determining the resonant coupling is formulated, 
that is consistent with the results of both the scattering method and the group analysis. In addition, 
a new type of photonic waveguide is proposed that manipulates the mechanism of photon hopping 
between photonic molecules and offers a new optical feature of twin waveguiding bandwidths. 

PACS numbers: 42.70.Qs, 42.82.Et, 42.60.Da, 71.15.Ap 



In the past decade, photonic defects have attracted 
much attention owing to their scientific and technological 
applications in the realization of high-Q microcavities or 
high transniittance waveguides (WGs) [1, 2, 3, 4, 5, 6, 7]. 
A defect atom can be embeded in a photonic crystal by 
perturbing the dielectricity of a selected crystal "atom" , 
that has photons with certain frequencies locally trapped 
within the band gap of the surrounding crystal struc- 
ture. If defect atoms are embeded by design to form the 
so called line defect then, within the band gap, it may 
provide a mechanism of photon propagation via hopping 
from one defect to its neighbors with a high transmit- 
tance [8, 9, 10]. Consequently, the integrated optical 
circuits of functional elements can be realized through 
skillful integration of photonic defects and WGs, and is 
expected to offer potential applications in telecommuni- 
cation [11, 12]. 

A point of importance but of much less attention is 
that the formation of photonic WGs are conventionally 
considered as arranging the desired defect atoms along a 
line, however, this approach has limited the potential of 
development. We suggest a broader vision of the manip- 
ulation of photonic defects through the investigation of 
photonic molecules that are defined as follows. In pho- 
tonic crystals, the defect atoms are closely arranged to 
form a structure that is similar to a real molecule — So 
named because, within the major band gap, the resonant 
modes confined by such a structure are closely analogous 
to the ground states of the molecular orbitals (MOs) of 
the real molecule. For example. Fig. 1 (a) shows a pho- 
tonic molecule named as the photonic benzene, whose 
defect modes are shown in Fig. 3. By employing the 
variational theory, the constraint determining the reso- 
nant coupling of photonic molecule is formulated, that is 
consistent with the results from those of both the scat- 
tering solution and the group analysis. In particular, 
manipulating the mechanism of photon hopping between 
photonic benzenes can provide the function of guiding 
photons along the benzene chain with a high transmit- 
tance, and presents a new optical feature of twin waveg- 
uiding bandwidths, as shown as Fig 4. 

Because of the importance on interpreting the mecha- 



nism of defect coupling, there are mainly two solid-state 
theoretical approaches, the tight-binding (TB) [8, 9, 13] 
and Wannier function methods [14, 15, 16], have been ap- 
plied to study the coupled cavities. The photonic version 
of TB method extends the idea of linear combination of 
atomic orbitals (LCAO), in which the defect modes are 
analogous to the atomic wave function, and suppose that 
only the nearest-neighbor couplings are relevant to find 
the dispersion relation for waveguide mode. For the lat- 
ter, the localized defect modes are expanded by Wannier 
functions to calculate their intensity variations, where the 
Wannier functions are calculated by plan- wave method or 
TB method coupled with supercell approximation. Here, 
another powerful approach of variational analysis is in- 
troduced for many defect-atoms system. 

Considering first an electric resonant mode Ed(r) of a 
single defect in a finite-size photonic crystal, the Maxwell 
equations obeyed by Ed(r) can be further simplified as 



HErf(r) =e<i(r)(c^<i/c)2Erf(r), 



(1) 



where the operator H is defined as — for the 2D system 
or V X Vx for the 3D system. Also, ed(r) denotes the 
dielectric constant of the single defect system, and LOd is 
the eigenfrequency of the eigenmode Erf(r). Those modes 
occure within the major band gap are most concerned in 
this paper, and can be taken as real and normalized by 

< Erf(r)led(r)lEd(r) / ed(r)Ed(r)-Ed(r)dr = 1. 

J domain 

For the same photonic crystal considered in (1) but 
embedded with a photonic molecule, the allowed resonant 
modes are assumed as a superposition of the individual 
defect mode. Basically, it is analogous to the idea of 
LCAO, namely 



nd 



E„(r) = _^C„,Ed(r-R,), 



(2) 



where E„(r) is the nth resonant mode of the photonic 
molecule, the number of defect atoms, C„i the unde- 
termined coefficient, and the coordinates of the ith 
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defect atom. Similarly, E„(r) satisfies (1) but with ed(r) 
replaced by the dielectric constant epm{j^) of the photonic 
molecule system, and replaced by the frequency a;„ of 
the eigenmode E„(r). That is 



HE„(r) = ep„(r)(w„/c)2E„(r). 



(3) 



Equation (2) associated with (3) is a linear variational 
problem. Assigning different coefficient C„i to each mode 
Ed may create different En, but the structure of photonic 
molecule will decide which resonant modes arc allowcxl. 
This inference will be reflected on restricting C„j to sat- 
isfy the minimum of functional frequency, defined as 



;E„(r) H E„(r) 



(E„(r)|ep„(r)|E„(r)) 



Namely, is equivalent to the familiar Rayleigh-Ritz 
principle. Here Hjj and Sjj denote, respectively, the ele- 
ments of the Ifamiltonian matrix and the overlap matrix. 
According to (1) and (2), can be written as 



Hy = (Ed (r - Ri) H Ed (r - R^) 



= < 



' {(^d/cf = a, for i = j; 

{uid/ cf f^i, for (i, j) being the nearest-neighbor; 



{uid/c) P2, for being the second-neighbor; 
^ (wd/c)^ /^s, for being the third-neighbor, 



where %) =< E<i(r - R,)|ed(r - Rj )|Ed(r - Rj) > de- 
notes the hopping integral, whose magnitude measures 
the coupling strength and decays rapidly with increas- 
ing the distance |Ri — Rj|, i.e. > l/^a] > l/^sl (the 
more the defect sites split, the weaker their coupling [17]). 
Therefore, hopping terms can be classified according to 
the separation of the coupled defects. Here, only three 
relevant hopping terms are considered. Moreover, un- 
der the assumption that each individual defect mode is 
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FIG. 1: Formation of (a) photonic benzene, (b) benzene WG 
of TT— type and (c) benzene WG of cr— type in a 2D hexagonal 
crystal; each solid dot denotes a defect atom of radius p. 



highly localized around its defect site, the field overlap 
between different defect atoms is small and the overlap 
integral can thus be approximated as 



Sij =< Ed(r - Ri)|ep„(r)|Ed(r - R^) >» Si, 



(4) 



Now, vary Cni fo minimize the functional frequency 
n, with the necessary condition of dCl/dCni = 0, i = 
1,2,..., Ud- One can obtain 



Cnj — 0. 



The constrain of resonant frequencies can thus be derived 
from the solvability condition of Cni- This leads to 



det 



Hy - {U>n/c) Si_ 



= 0. 



(5) 



Equation (5) indicates that the allowed resonant fre- 
quencies in a given photonic molecule are dominated by 
hopping integrals. Furthermore, these hopping terms are 
dependent upon the dielectric structure of the photonic 
molecule. Therefore, every resonant mode is character- 
ized by photonic molecule and exhibits different optical 
transmittance. To check the accuracy of Eq. (5), we first 
apply Eq. (5) to the photonic benzene, that yields 



0, (6) 



where we let 7 = (a;„/c)^ and /3j = {cOd/c)'^(3i for sim- 
plification. The determinant in Eq. (6) is called the 
6th-order circulant, and is equivalent to 
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0, 



where ei,e2...e6 are the six roots of unity, i.e. e„ = 
exp{2TTin/Q). Hence, for n = 1 . . .6, the frequencies a;„ 
of the six resonant modes E„ can be found as 

: Wd\/ 1 + 2 cos — /3i + 2 cos — /?2 + (-l)"/33- 

(7) 

Obviously, Eq. (7) indicates that there are two doubly- 
degenerate frequencies with n = 1,5 and n = 2,4, and 
two nondegenerate frequencies with n = 3 and n = 6, 
thus, four high-Q resonant frequencies will occur within 
the major band gap. From the view point of the sym- 
metry group , the photonic benzene belongs to the point 
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group Cgt, for the 2D or Dgh for the 3D systems, whose 
irreducible representation F on defect sites can be re- 
duced to the decomposition T = Ai + B2 + Ei + E2 or 
r = + + Eig + E2U, respectively. Exactly speak- 
ing, it again illustrates two doubly degenerate modes of 
El (or Eig) and E2 (or E2u), and two nondegenerate 
modes of Ai (or A2u) and B2 (or B2g)- Furthermore, 
these predictions are also consistent with the numeri- 
cal solution of Eq. (3), that is calculated by scattering 
method (cf. [17]). The resultant transmission for the 
2D case is plotted in Fig. 2, and the allowed resonant 
modes are shown in Fig. 3, where we consider a 2D finite- 
size hexagonal photonic crystal with a dielectric contrast 
ratio 8.4/1.0 (rod/background) and a radius-to-spacing 
ratio 0.4/4.0. The defect radius p is taken as zero. 

Figure 2 shows four nomalizcd resonant frequencies 
0.369, 0.375, 0.455 and 0.508, that split from cOd = 0.419 
due to the defect coupling in a photonic benzene. By sub- 
stituting these five frequencies into Eq. (7), the hopping 
terms (32 and /Js can be calculated by least-squares 
method, and they are 0.178, 0.051 and -0.010. Table 
I summariz(^s four cases of defect sizes, in which the all 
values of hopping parameters achieve the accuracy of two 
decimal places for Eq. (7). It clearly shows that the 
larger the defect radius is, the smaller hopping parame- 
ters. This means that the defect couplings in photonic 
benzene become weaker as the defect radius is increased. 
Moreover, when the defect radius p is increased up to 
about 0.2, the property of four transmission peaks disap- 
peares, owing to that the shallow perturbation of dielec- 
tricity for defect atoms will create shallow modes [1] . 

In fact, Eq. (3) is equivalent to the effective Shrodinger 
equation of Hitckel tt— electron theory (developed in 1931 
[18]), if we make the resonant modes of 3D photonic 
benzene be analogous to the tt— electrons of benzene 
molecule. However, Eq. (3) is much simpler. Theoret- 
ically, the TT-electrons arise from a planar unsaturated 
organic molecule whose MOs can be divided into the 
a and tt MOs according to the reflection symmetry in 
the molecular plane. Both systems belong to the same 
point group De/i and have the same degeneracy (cf. [19], 
p. 261), but possess completely different meanings. Con- 
ceptually, the photonic molecule acts as a perfect model 
of artificial molecule, since the resonant modes are much 
easier to be realized than the bonding orbitals of real 



TABLE I: The calculated values of normalized resonant fre- 
quencies and hopping parameters for photonic benzene with 
different 

p LUd CJl UI2 L^-A 1^4 01 02 03 

0.0 0.419 0.455 0.375 0.369 0.508 0.178 0.051 -0.010 

0.05 0.417 0.452 0.373 0.367 0.505 0.177 0.051 -0.009 

0.1 0.409 0.442 0.369 0.361 0.493 0.171 0.048 -0.006 

0.15 0.395 0.423 0.361 0.351 0.468 0.155 0.040 -0.002 
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FIG. 2: Transmission amplitude of electric field as function 
of normalized frequency for a photonic benzene (solid curve) 
and a defect atom (dotted curve). The marked values are 
four frequencies of the high-Q modes which are split from the 

ujd due to defect coupling and fallen in the major band gap 
whose rang is as shown in Fig 4. 



molecule. Similar phenomena can also be found in the 
quantum-dot molecules [20] or the coupled pairs of GaAs 
cavities (note that these systems are also termed as pho- 
tonic molecule, cf. Bayer et al. [21]). Fig. 3 shows the 
lowest resonant modes allowed in a 2D photonic benzene 
with p = 0.0, and they arc labelled according to the sym- 
metry species of the group C^y, which are analogous to 
the six TT MOs of benzene molecule but is lacking of the 
Cih symmetry. 

Most importantly, by applying the modular concept 
of photonic benzene can create photonic WGs, and we 
called these WGs as benzene WGs. In the chemistry 
terminology, benzene WGs can be classified as tt and a 
types corresponding to the bonding types between two 
real benzene rings, as illustrated in Fig. 1 (b) & (c). 
It is remarkable that transmission of the benzene WGs 
reveals the special feature of a twin waveguiding band- 
widths marked as Auja and Acufj, where Aw^ < 
as shown in Fig. 4. This means that the benzene WGs 
are able to provide two working bandwidths at the same 
time. In addition. Fig. 5 shows that a TM light with the 
E2 mode travels through a 90" bend from the tt type to 
the a type. Of course, the same phenomenon can also be 
observed in other modes. 

In conclusion, we siiggcst a new and practicable idea 
for the manipulation of photonic defects, that includes 
the so-called photonic molecule and benzene WG. The 
optical properties of photonic molecules has been inves- 
tigated by variational theory, which shows that the al- 
lowed resonant frequencies inside a photonic molecule 
are dominated by hopping parameters through the con- 
straint (5). Taking the photonic benzene as an example, 
six resonant modes with two doubly-degenerate and two 
nondegenerate arc found and verified by both of the scat- 
tering method and group theory. Especially, the benzene 
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WGs created by the modular manipulation of photonic 
benzenes are demonstrated to possess the interesting fea- 
ture of a twin waveguiding band widths. Namely, benzene 
WGs provide two working bandwidths at the same time 
and make the function of guiding photons more flexible. 

This work is supported in part by the Global Fiberop- 
tics, Inc. 
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FIG. 4: Transmission created by benzene WGs of (a) tt— type 
and (b) a— type with p = 0.0. Both types present the feature 
of twin waveguiding bandwiths, marked as Auja and Aufj. 
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FIG. 3; Resonant electric field pattern in a photonic benzene 
of p = for the six lowest resonant modes with normalized 
frequencies of (a) 0.369 with Ei mode, (b) 0.369 with Ei 
mode, (c) 0.455 with E2 mode, (d) 0.455 with E2 mode, (e) 
0.375 with Ai mode and (f) 0.508 with B2 mode, respectively. 
In general, the more node, the higher frequency. 



FIG. 5; Eletric field distribution of a TM light with E2 mode 
travels through a 90° bend in a benzene WG where the junc- 
tion connects a tt type WG with a a type; incident wave with 
' normalized frequencies 0.449 enters from the upper left. 



B 61, R11855 (2000). 
[11] H. Kosaka, T. Kawashima, A. Tomita, M. Notomi, T. 

Tamamura, T. Sato, and S. Kawakami, Appl. Phys. Lett. 

74, 1370 (1999). 
[12] A. R. McGurn, Phys. Rev. E 65, 075406 (2002). 
[13] E. Lidorikis, M. M. Sigalas, E. N. Economou, and C. M. 

Soukoulis, Phys. Rev. Lett. 81, 1405 (1998). 
[14] K. M. Leung, J. Opt. Soc. Am. B 10, 303 (1993). 
[15] J. P. Albert, C. Jouanin, D. Cassagne, and D. Bertho, 

Phys. Rev. B 61, 4381 (2000). 
[16] A. Garcia-Martin, D. Hermann, K. Busch, and P. Wolfie, 

Mat. Res. Soc. Symp. Proc. 722, L 1.1.1 (2002). 
[17] G. Tayeb and D. Maystre, J. Opt. Soc. Am. A 14, 3323 

(1997). 

[18] E. Hiickel, Z. physik 70, 204 (1931). 

[19] P. W. Atkins and R. S. Friedman, Molecular Quantum 

Mechanics (Oxford University, Oxford, 1997). 
[20] M. Pi, A. Emperador, M. Barranco, F. Garcias, K. Mu- 



5 



raki, S. Tarucha, and D. G. Austing, Phys. Rev. Lett. L. Rcincckc, P. A. Knipp, A. A. Drcmin and V. D. Ku- 

87, 066801 (2001). lakovskii, Phys. Rev. Lett. 81, 2582 (1998). 

[21] M. Bayer, T. Gutbrod, J. P. Reithmaier, A. Forchel, T. 



